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The effect of dc electrical stress and breakdown on Josephson and quasiparticle tunneling in 
Nb/Al/AlO^/Nb junctions with ultrathin AlOa; barriers typical for applications in superconductor 
digital electronics has been investigated. The junctions' conductance at room temperature and 
current-voltage (I-V) characteristics at 4.2 K have been measured after the consecutive stressing 
of the tunnel barrier at room temperature. Electrical stress was applied using current ramps with 
increasing amplitude ranging from to ~ 10007 c corresponding to voltages across the barrier up to 
~ 0.65 V, where I c is the Josephson critical current. A very soft breakdown has been observed with 
polarity-dependent breakdown current (voltage). As the stressing progresses, a dramatic increase in 
subgap conductance of the junctions, the appearance of subharmonic current steps, and gradual in- 
crease in both the critical and the excess currents as well as a decrease in the normal-state resistance 
have been observed. The observed changes in superconducting tunneling suggest a model in which 
a progressively increasing number of defects and associated additional conduction channels (super- 
conducting quantum point contacts (SQPCs)) are induced by electric field in the tunnel barrier. 
By comparing the I-V characteristics of these conduction channels with the nonstationary theory 
of current transport in SQPCs based on multiple Andreev reflections by Averin and Bardas, the 
typical transparency D of the induced SQPCs was estimated as D ~ 0.7. The number of induced 
SQPCs was found to grow with voltage across the barrier as sinh(VyVb) with Vo = 0.045 V, in 
good agreement with the proposed model of defect formation by ion electromigration. The observed 
polarity dependence of the breakdown current (voltage) is also consistent with the model. Based on 
the observed magnitude of breakdown currents, electric breakdown of Aldbarrier during plasma 
processing was considered to be an unlikely cause of fabrication-induced, circuit pattern-dependent 
nonuniformities of Josephson junctions' critical currents in superconductor integrated circuits. 



I. INTRODUCTION 

Aluminum oxide (AlO^) is widely used as a barrier ma- 
terial in various applications involving tunnel junctions. 
It was also considered as a potential gate oxide in ad- 
vanced memory devices and metal-oxide-semiconductor 
(MOS) transistors. Dielectric reliability issues such as ox- 
ide barrier stability, leakage currents, and electric break- 
down are very important for electronics applications, es- 
pecially for magnetic tunnel junctions (MTJs) used for 
magnetic random access memories and superconducting 
tunnel junctions (STJs) used for superconductor digital 
circuits requiring ultrathin (~ 1 nm) tunnel barriers. The 
physics of dielectric breakdown in ultrathin barriers is of 
great interest in its own right. Oxide breakdowns are 
usually classified into two modes: intrinsic and extrin- 
sic. Though the difference is somewhat blurry, extrinsic 
breakdowns are those caused by defects introduced or 
created during oxide growth, whereas intrinsic ones are 
the property of a perfect dielectric. In relatively thick 
oxide layers such as those used as gate dielectric in MOS 
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transistors, electric breakdown usually proceeds by the 
accumulation of defects (traps) inside the dielectric until 
a percolation pass is formed, at which point the resis- 
tivity suddenly decreases from a very high value to a 
very low value, and a hard breakdown occurs. [J, [|J In 
thinner oxide layers the breakdown often has a soft char- 
acter which is characterized by small gradual changes in 
resistance. [3, E| The thickness of ultrathin tunnel barri- 
ers used in superconductor electronics is a couple of oxide 
monolayers, only a few interatomic distances. Therefore, 
any defects formed in the oxide as a result of electrical 
stress (e.g., displaced ions, oxygen vacancies, etc.) should 
create additional conduction channels with significantly 
increased transmission probability and consequently dra- 
matically alter the quasiparticle and Cooper-pair tunnel- 
ing. In other words, the percolation path forming at 
breakdown may consist of just a single defect (trap) and 
hence the breakdown may be very soft. 

There have been several publications on the reliabil- 
ity and breakdown of aluminum oxide layers with thick- 
nesses above ~ 3 nm as a new gate oxide and ~ 1 nm in 
MTJs.fl, IE 0, H, Q The existence of both intrinsic and 
extrinsic breakdown modes was suggested. Q The intrin- 
sic mode was associated with a hard breakdown, and was 
suggested to be related to the chemical bond breaking in 
applied electric field. [IE HH, Q3 The extrinsic mode was 
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associated with a soft, gradual breakdown. It was sug- 
gested to be related to pre-existing pinholes in the barrier 
which grow in area as breakdown progresses due to Joule 
heating and/or electric field effect. 

Superconductor-insulator-superconductor (SIS) junc- 
tions offer unique opportunities in studying breakdown 
mechanisms in ultrathin oxides because both the quasi- 
particle and Josephson tunneling in STJs are extremely 
sensitive to the barrier properties and boundary condi- 
tions at the metal-oxide interfaces. In contrast to MTJs, 
pre-existing pinholes in SIS junctions are easily iden- 
tifiable because they carry supercurrent thus creating 
nonuniform Josephson current distribution and dramat- 
ically increasing subgap conductance. Whereas a mi- 
croshort in MTJs was associated with a junction hav- 
ing the resistance-area product RA of ^0.8 ft ■ /mi 2 , 
and no tunneling magnetoresistance,Q STJs with even 
lower values of RA product exist and display interesting 
Josephson tunneling properties. [IH 0, El 

Superconductor digital electronics utilizing SIS junc- 
tions has a potential for sub-THz clock frequencies and 
ultra-low power dissipation for digital signal processing, 
high-performance communications and computing, [la] 
Recently, complex superconducting circuits based on 
Rapid Single Flux Quantum (RSFQ) logic such as 
analog-to-digital converters and digital RF receivers con- 
taining thousands of logic gates with clock frequen- 
cies ^30 GHz have been demonstrated, operating not 
only in liquid He but also on commercial closed-cycle 
crvocoolers. (l7j Increasing the clock frequencies of super- 
conductor integrated circuits to ~ 100 GHz would require 
employing high- J c junctions with RA products below ~1 
£1 ■ /im 2 , perha ps t he thinnest tunnel barriers among all 
known devices.]!!, 

Dielectric reliability may not appear to be important 
for superconducting digital circuits because they oper- 
ate at very low temperatures and at very low voltages 
(~1 mV). Its significance however arises from the possi- 
bility that tunnel barrier degradation may occur during 
integrated circuit fabrication. For instance, the current 
state of the art in RSFQ circuits has been plagued by 
limited circuit yield brought about to a large extent by 
fabrication-induced variations on the J osep hson critical 
current (7 C ) of the tunnel junctions. [H, These varia- 
tions may be related to dielectric barrier degradation and 
electrical breakdown. It has been observed, for example, 
that the J c of Josephson tunnel junctions may depend on 
how the junction is wired to other circuit elements, and 
in particular, at which step in the fabrication process 
the junction makes electrical contact with the circuit's 
ground plane. [2CJ For series arrays of nominally identi- 
cal tunnel junctions, it was frequently observed that the 
Josephson critical current of a few junctions (usually of 
the first and the last junction in the array) is significantly 
larger than for the rest of the junctions. This cannot be 
simply explained by a variation in the area of that one 
junction coming from the lithography and etch processes 
of junction definition. Neither can it be explained by 



a random fluctuation in the tunnel barrier transparency 
in that particular junction, considering that the effect re- 
produces in different arrays and the junctions in the array 
are just ~10 fim apart. Instead, it was suggested that 
the above phenomena are brought about by electrical 
currents flowing through the tunnel barriers, a result of 
plasma processing steps which follow the SIS trilayer de- 
position. Recent experiments involving tunnel junctions 
protected from plasma process-induced electric stress by 
current-limiting resistors support this suggestion. [2lj 

Surprisingly, there has been almost no research on 
the reliability of ultrathin AlO^ barriers in STJs, except 
for early works on Al/AlO^/Pb junctions which studied 
the effects of electric annealing on the barrier thickness, 
height and asymmetry of AlO^ barriers formed by plasma 
oxidation. [22j These junctions, however, have no practi- 
cal application, and the changes in superconducting and 
Josephson properties were not studied. 

In this paper the effect of applied dc electrical 
stress on quasiparticle and Josephson tunneling in 
Nb/Al/AlOz/Nb junctions was investigated. The study 
focused on this type of junctions because of their domi- 
nant use in superconductor digital and analog electronics. 



II. FABRICATION 

The Nb/Al/AlO^/Nb junctions used in this study were 
fabricated at HYP RES, Inc. using an 11-level process for 
superconductor integrated circuits. [H, [23[ The fabrica- 
tion was performed on 150-mm Si wafers. The process is 
based on in-situ Nb/Al/AlO^/Nb trilayer deposition. [2~j] 
Specifically, Nb/Al bilayer (150 nm and 8 nm, respec- 
tively) deposition is followed by AlOa; formation by room 
temperature oxidation of Al. The AlO^ layer is then 
topped off by the deposition of the Nb counter-electrode 
(50 nm) . All metal layers are deposited by dc magnetron 
sputtering in a cryopumped vacuum system with base 
pressure of 1 x 10~ 9 Torr. Two Josephson critical cur- 
rent densities, J c (1 kA/cm 2 and 4.5 kA/cm 2 ), were tar- 
geted, obtained by Al oxidation for 15 min at oxygen 
pressure of 170 mTorr and 18 mTorr, respectively. After 
the counter-electrode etch process that defines the junc- 
tions, their interior was sealed along the perimeter and 
sidewalls by an anodization layer composed of AI2O3 and 
Nb205 in order to protect the barrier from reacting with 
process chemicals used in subsequent fabrication steps. 

Circular JJs with design radii of 2.00 /im (A = 12.6 
/im 2 ) and 0.95 fim (A = 2.8 /im 2 ) for wafers with J c = 1 
kA/cm 2 and 4.5 kA/cm 2 were arranged on 5 x 5-mm 2 
chips referenced according to wafer number and the co- 
ordinates (in units of 5 mm) of the chip location on the 
wafer. For example, a chip from wafer KL1004 with lo- 
cation (-5,7) is to be called KL1004N5P7. 
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III. EXPERIMENT 

Each chip was mounted inside a magnetically shielded 
cryoprobe. The junctions were measured using a low-pass 
filtered four-probe setup with a Keithley 2000 voltmeter 
and a Keithley 6220 current source. Electrical stress was 
applied at room temperature by ramping current through 
the junction up to a preselected value Is and back down 
to zero. The stress application was preceded and followed 
by a measurement of the junction's room-temperature re- 
sistance, using a low current of 500 /iA. The tunneling 
I-V characteristics were then measured in the supercon- 
ducting state at T — 4.2 K with the cryoprobe submerged 
in a liquid He. The next stress/measurement cycles were 
then performed using progressively higher values of Is ■ 

The effect of stress polarity was also investigated. Here 
we define positive stress as current flowing from the 
counter electrode (positive potential on the top Nb layer) 
to the base electrode (Nb/Al bilayer), while negative 
stress is current flowing in the opposite direction. At 
room temperature, there is a resistance in series with the 
tunnel barrier associated with the normal resistance of 
interconnects to the junction. This series resistance was 
estimated from the interconnects' geometry using a sep- 
arately measured sheet resistance of the layers involved. 
This resistance was assumed to remain unchanged by 
stress applications and was simply subtracted from the 
total measured resistance, thus allowing us to estimate 
the potential difference that develops across the tunnel 
barrier during stressing. Parameters of some of the stud- 
ied junctions are given in Table HI 

IV. EXPERIMENTAL RESULTS 

Fig. Q] shows the resistance at low currents (< 500 
/j,A) of four Nb/Al/AlO^/Nb junctions after each subse- 
quent stress application. The junction resistance remains 
roughly constant after stressing with low currents until 
a threshold stress current for breakdown Isb is reached, 
above which the resistance starts to decrease indicating 
irreversible changes in the barrier. Each succeeding de- 
crease in resistance is apparently a cumulative effect of all 
previous stress applications. The threshold current varies 
among nominally identical junctions on the wafer and 
from wafer to wafer, suggesting some statistical nature of 
the barrier breakdown. Despite these variations, the ob- 
served threshold current for the positive stress (current 
from Nb counter-electrode to Al) is consistently higher 
than for negative stress as listed in Table fl] If the re- 
sistance after stressing is scaled with the resistance in 
the initial junction and the stress current is scaled with 
the threshold current, all the curves in Fig. [1] collapse 
onto a single curve (see Fig. [1] inset) suggesting a uni- 
versal breakdown mechanism and a universal character 
of resistance changes due to electrical stress. 

The irreversible decrease of the tunnel barrier resis- 
tance may be due to (a) a decrease in the average bar- 
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FIG. 1: Junction resistance at room temperature after elec- 
trical stressing of the junctions in Table [I] (+)/(-) indicates 
positive/negative stress polarity. Despite small differences be- 
tween individual junctions, the stress-induced irreversible re- 
sistance changes are very similar as demonstrated in the Inset 
showing the normalized resistance R/R n at room temperature 
as a function of normalized stress current Is /Isb for the same 
junctions (see Table U]). 



rier height and/or thickness; (b) the barrier becoming 
nonuniform due to formation of additional conduction 
channels (regions with increased barrier transparency 
which are often called micro- or nano-shorts); and (c) 
a combination of the above. However, room tempera- 
ture measurements alone are insufficient to distinguish 
between these possibilities and measurements in the su- 
perconducting state of the junction electrodes are needed. 

The electric stress-induced changes in the tunnel bar- 
rier properties are clearly seen in the Josephson I-V char- 
acteristics at T — 4.2K shown in Fig. EJa). Five main 
features are worth mentioning. First, the Josephson crit- 
ical current, I c (defined here as the switching current) 
increases with the stress current. Second, the normal- 
state resistance, R n of the junction decreases with the 
stress current. Third, the so-called knee structure in the 
quasiparticle tunneling just above the gap voltage of the 
junction also shifts to higher currents as the stress in- 
creases. Fourth, the junction conductance in the sub- 
gap region of voltages increases with the stress current. 
The retrapping current also increases with the increase 
in subgap conductance. Finally, the gap voltage V g is 
nearly independent of the stress current and decreases 
only slightly at very high stress currents. The increase in 
the subgap conductance is by far the most pronounced 
change in the I-V curves of the stressed junctions. For 
instance, after applying a stress current of 96 mA, the 
Josephson critical current of the junction increases by a 
factor of ~ 2.5 whereas the subgap conductance at 2 mV 
increases by a factor of ~ 25. 
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TABLE I: Summary of initial resistances, initial critical current, and breakdown stress current for the samples shown in Fig. [T] 
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!: (a) I-V characteristics of Nb/Al/AlO^/Nb junction 
4.2 K (initial J c = 1 kA/cm 2 ) after each application of 
electrical stress show increasing 7 C , decreasing R n , increasing 
subgap conductance, and increasing excess current. The gap 
voltage Vg = 2A/e and the current step at V = V s remain 
almost unaffected by electric stress in the wide range of stress 
currents from up to ~ 27sb. At higher stress currents the 
gap structure broadens and diminishes, and at Is ~ 37sb the 
junction loses all remaining signatures of the tunnel junction. 
Numbers in the legend indicate the applied positive stress 
current in mA and identify the curves from top to bottom, (b) 
Blow-up of the return branches of I-V curves of KL1004N5N6 
after each stress application clearly shows the development of 
current steps (subgap structure) at subharmonics of the gap 
voltage. 
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FIG. 3: I c Rn product (left scale, dotted lines) and i? ag at 2 
mV (right scale, solid lines) values after each stress applica- 
tion normalized to their initial values I c oRnO and 7? sg o in the 
unstressed junctions. Parameters of the junctions are given 
in Table |U 



V. DISCUSSION 

The experimental results presented in the previous sec- 
tion show that there is a broad range of electric stress 
currents through Nb/Al/AlO^/Nb junctions (and corre- 
sponding voltages across the AlO^ tunnel barrier) where 
a gradual change of the barrier properties occurs. This 
can be interpreted as a soft breakdown of the oxide bar- 
rier in contrast to the sudden changes (hard breakdown) 
usually observed in thicker oxide layers of SiC>2 and AlO^ . 
For the junctions with J c = 1 kA/cm 2 , this range spans 
from the typical threshold current I§b ~ 30 mA (cor- 
responding to a stress current density Js ~ 2.4 x 10 5 
A/cm 2 ) to I s ~ 120 mA ( J s ~ 10 6 A/cm 2 ) above which 
the stressed devices lose all the signatures of supercon- 
ducting tunnel junctions. For the positive stress polar- 
ity, the threshold current in Table HI Isb = 56 mA cor- 
responds to the voltage across the barrier = 0.519 
V. For the negative stress polarity, Isb = 32 mA corre- 
sponds to VC = 0.252 V. The typical barrier thickness d 
for the Josephson current densities used in this work is 
~ 1 nm or less.[25|, [2(| Hence, the typical electric fields 
across the AlO^ barrier at which the irreversible changes 
start are = 2.52 x 10 s V/m and E+ = 5.19 x 10 8 
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FIG. 4: A comparison of the I-V curve of a stressed junction 
KL1023N5P8 (dashed curve) with post-stress Josephson crit- 
ical current density J c ~ 9 kA/cm 2 (initial J c = 4.5 kA/cm 2 ) 
and an as-fabricated, unstressed junction KL1013N1P1 (solid 
curve) with J c ~ 11 kA/cm 2 . Although the as- fabricated 
junction has an even larger J c (larger average barrier trans- 
parency), its I-V curve is very different from the electrically 
stressed junction: it has no appreciable subgap conductance, 
no subharmonic current steps, and no excess current; its 7 c 7? n 
product and the current step at V g are close to the values 
given by the microscopic theory for tunnel junctions with low 
transparency. Presumably, the as-fabricated junctions have 
a uniform tunnel barrier whereas the barrier in electrically 
stressed junctions becomes nonuniform. The I-V curve of the 
initial, unstressed junction KL1023N5P8 is also shown (dot- 
ted curve). 




1 2 

Voltage (mV) 

FIG. 5: The return branch of the I-V characteristic (dotted 
curve) of KL1004N5N6 after application of 80 mA stress along 
with differential conductance dl/dV (solid curve) showing 
peaks corresponding to multiple Andreev reflections (MAR) 
of quasiparticles. Inset shows the voltages corresponding to 
conductance peaks; the straight line is a fit to 2 A/en (n = 
1,2,.. .) dependence expected for MAR, giving 2A/e = 2.696 
mV. 



V/m for the negative and positive stress polarities, re- 
spectively. 

The difference between V^~ and VC~ was also observed 
in AlOa; breakdown measurements in magnetic tunnel 
junctions[l, [j| where it was speculated to be a result 
of different surface roughness of the two metal/oxide in- 
terfaces. We suggest that the difference between V^~ 
and V^~ is more fundamental and is a result of differ- 
ence in work functions of the junction electrodes. Indeed, 
from the contact potential difference measurements, the 
work functions in the electrodes are (fA\ — 4.19 eV and 
VNb = 4.37 eV for Al and Nb, respectively. Therefore, 
in a Nb/Al/AlO^/Nb junction the potential barrier be- 
comes asymmetric and there is an internal electric field 
across the barrier E ult = A</? Nb /Ai/ed ~1.8 x 10 8 V/m 
directed from the Nb/Al base electrode towards the Nb 
counter electrode and arising from the difference in work 
function A<^ N b/AZ — fNb — <PAl — 0.18 eV. This bar- 
rier asymmetry agrees well with the result obtained from 
the asymmetry of conductance vs. voltage characteris- 
tics of Nb/Al/AlOz/Nb junctions with low tunnel bar- 
rier transparency. pHi l26j For the negative stress polarity 
(higher potential on Nb counter-electrode), the external 



field adds up to the internal field E t 



E,-, 



the same direction, whereas for the positive stress po- 
larity the external field is opposite to the internal one 
£ t + = E cx t — i?i n t, where E t is the total field in the di- 
electric and E cyA is the externally applied electric field. 
If defects in the barrier start to form when the net inter- 
nal field reaches some critical level E c (e.g., as a result of 
ion electromigration) , from the onset of irreversible resis- 
tance changes (soft breakdown), we estimate the critical 
voltage to be V c = (V b + + V b _ )/2 = 0.385 V, correspond- 
ing to E c w 3.85 x 10 s V/m. The barrier asymmetry 
in this model corresponds to — V b ~)/2 = 0.13 V, in 
excellent agreement with the 0.18 V difference between 
work functions. 

Let us see what can be inferred about the nature of 
stress-induced changes to the barrier from the data pre- 
sented in Sec. IIVI Clearly the increase of the critical 
current and decrease in the normal resistance in stressed 
junctions indicate that the critical current density and 
the average barrier transmission increase. The question 
is whether the tunnel barrier remains uniform or not. If 
the barrier remains uniform and only its height and/or 
thickness is changing with the applied stress, one would 
expect the junction I c R n product to remain constant be- 
cause in the limit of low barrier transmission probability, 
it is given by the Ambegaokar-Baratoff (AB) relationship 
IcRn — 7rA/2e for T <C X^,]^ and the gap voltage was 
found to remain nearly constant after barrier stressing. 
Contrary to this, the measured I c R n product increases 
significantly in the stressed junctions as can be seen in 
Fig. [31 This may indicate a transition from tunnel junc- 
tion to point contact behavior because for point contacts 
between two superconductors at T <C T c the I c R n prod- 
uct is greater than in the tunnel junction by a factor 
of 1.32 in the dirty limit and by a factor of 2 in the 



6 



clean limit, as was shown by Kulik and Omel'yanchuk 
(KO).[H, III It was also found by Arnold [Ij that, 
for an insulating and structureless barrier of arbitrary 
thickness, the I c R n product in general depends on the 
barrier transmission probability D, and reduces to the 
AB result at D -C 1 and to the KO result for the clean 
point contact at D = 1. 

It is easy to show, however, that the barrier does not re- 
main uniform and structureless after stressing. For this, 
we compare the I-V curves of a stressed junction and an 
as-prepared, unstressed junction with the same average 
Josephson critical current density. The as-prepared junc- 
tion is presumably as uniform as possible for the given 
fabrication method. The I-V curves are shown in Fig. [4] 
The stressed junction with J c ~ 9.3 kA/cm 2 clearly has 
different I-V characteristics (especially at V < V s ) than 
the as-prepared junction with an even higher J c ~ 11 
kA/cm 2 . This comparison indicates that the changes 
in the I-V curves and the increase in the I c R n prod- 
uct brought about by electric stress cannot be explained 
by a gradual increase in transparency of a uniform and 
structureless barrier. A dramatic increase in the subgap 
conduction in the stressed junctions, the appearance of 
pronounced features at subharmonics of the gap voltage, 
2 A/ en at n — 2, 3, . . . shown in Fig. [5] along with the in- 
crease in / c 7? n strongly indicate that additional conduc- 
tion channels with increased transparency are gradually 
formed as a result of electric stress. These conduction 
channels can be viewed as point contacts between the 
junction electrodes, pinholes, or nanoshorts in the bar- 
rier, contributing to both the quasiparticle and Cooper- 
pair transport in parallel to the main tunnel barrier which 
remains largely unmodified by the stress. 

Subgap features similar to those appearing in electri- 
cally stressed junctions (Fig. [^b) and Fig. [5]) have long 
been observed in Nb/Al/AlO^/Nb and other types of 
junctions [H, Q [l5| with high critical current densi- 
ties and attributed to MAR. Two models were proposed 
to explain the appearance of MAR in high- J c junctions. 
Kleinsasser et al. [32] suggested that, as oxygen expo- 
sure during Al oxidation decreases, the barrier becomes 
nonuniform, consisting of regions with low transparency 
(good tunnel barrier) and regions with high transparency 
(pin holes). That is, the transparency distribution has 
two sharp peaks, one at a low D value and another one 
at D ~ 1. As oxygen exposure decreases, the relative con- 
tribution of the second peak increases, and so does the J c . 
Naveh et al. [33[ argued that the defects in the ultrathin 
tunnel barrier are naturally occurring and, therefore, the 
transparency distribution in high- J c junctions is univer- 
sal and given by the Schcp-Bauer (SB) distribution [33j 



p(D) = 



G 



1 



(1) 



They found that SB distribution gives a better fit to ex- 
perimental I-V curves than averaging over the Dorokhov 
distribution describing the transparency distribution in 
long disordered conductors — not a surprise though be- 



cause the Dorokhov distribution should not be applicable 
to short channels by definition. There is ample experi- 
mental evidence however that the behavior of high-J c 
junctions is not universal. For instance, the value of J c 
at which MAR steps become visible strongly depends on 
the junction fabrication procedure and varies between 
different experimental groups. [Til Il5l] Our experimental 
results show that MAR steps can be created at will in 
tunnel barriers with any initial transparency as a result 
of electrical stress. Therefore, what were perceived as 
natural and universal defects appearing in the oxide bar- 
rier as a result of a short oxidation process UK could 
be simply a result of damage to the very thin tunnel 
barrier induced during the junction (wafer) fabrication 
processes. This could explain why improvements in the 
junction fabrication process usually shift the J c level at 
which MAR steps begin to appear towards higher values. 

Although no microstructural characterization of the 
barrier changes was undertaken in this work, a plausi- 
ble scenario of how the additional conduction channels 
are induced by electric stress can be proposed. It is 
known that AlOa; formed by room temperature oxida- 
tion is amorphous and nonstoichiometric, i.e. contains 
a large amount of oxygen vacancies. [35| At the oxida- 
tion conditions used in our work, the barrier thickness 
is only ~ 1 nm, i.e. just about two to three nearest- 
neighbor distances. Therefore, a displacement of even a 
single atom from the barrier can significantly increase lo- 
cal transparency. If a dc electric field is applied to the 
oxide barrier, electromigration of cations (Al 3+ ) and an- 
ions (0 2 ~, — ) in opposite directions may occur. From 
anodic oxidation of Al in electrolytes it is known that the 
anodizing ratio for Al is ar — 1.3 nm/V, i.e., 1.3 nm of 
AlO x is formed for each volt applied to the Al/electrolyte 
interface. That is, on average, oxygen ions are trans- 
ported inside the aluminum electrode by ~ 1 nm per 
1 volt of applied voltage. Ionic current density in the 
oxide is a strong function of the electric field strength 
Ji = aexp((3E — 7-E 2 ).[36|, H3] In our case, there is no 
electrolyte providing a constant supply of oxygen atoms 
for the oxide to grow. However, some oxygen ions from 
A\O x barrier will electromigrate into Al underlayer of the 
Nb/Al/AlOz structure if the Nb/Al base electrode is at a 
positive potential (negative stress polarity in our nomen- 
clature) while Al 3+ cations should migrate towards the 
interface with Nb counter electrode. At a positive po- 
tential on Nb counter electrode (positive stress), oxygen 
ions from the AlO^ barrier will move into Nb counter 
electrode, depleting the barrier. The average oxygen dis- 
placement at the onset of breakdown can be estimated 
using the typical threshold voltage V~ — 0.25 V and 
ar = 1.3 nm/V, giving ~ 0.3 nm. This closely corre- 
sponds with 0-0 nearest-neighbor distance of 2.8 A in 
amorphous alumina. [38] 

Any ion displacement in the barrier under the effect 
of an applied electric field creates a defect and alters the 
barrier transparency. The length of the newly created 
conduction channel is about the same as the barrier thick- 
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Let us see what properties of conduction channels cre- 
ated in the AlO^barrier by electrical stressing can be in- 
ferred from the measured I-V curves shown in Fig. [5] 
Because of space constraints we will restrict our analysis 
to the positive stress polarity only. The case of the neg- 
ative stress polarity is completely analogous and will be 
presented elsewhere. First, we note that the amplitude of 
the current step A/ ss at V g = 2A/e practically does not 
change in stressed junctions although the subgap con- 
ductance dramatically increases. (In the first approxi- 
mation we will neglect some broadening of the current 
step at Vg and a small decrease in V g at extremely high 
stress currents which could result from the smearing of 
the metal- insulator interfaces.) In the microscopic theory 
of tunneling in superconductors, the size of the current 
step is proportional to the area of the tunnel barrier and 
is given by 



FIG. 6: Normalized I-V characteristics of the additional con- 
duction channels (SQPCs) created by positive electric stress- 
ing at currents 7s = 56, 64, . . . , 96 mA. The curves were ob- 
tained by subtracting the I-V characteristics of the initial, 
unstressed junction Io(V) from the I-V curve after each stress 
application, 7 S qpc(V) = 7(V) - Io(V). The dip at eV/A = 2 
in the curves is an artifact of the subtraction procedure due 
to some broadening and slight decrease of the gap after elec- 
tric stressing. Theoretical I-V curves for a single SQPC with 
varying transparency D = 0.01, 0.1, 0.2, 0.3, . . . , 0.8, 0.9, 0.99 
(bottom to top) from [i^] are also shown. As can be seen, all 
the obtained experimental dependences fall within the range 
of theoretical curves corresponding to 0.6 < D < 0.8. 



ness which is much less than the coherence length £ and 
the magnetic field penetration depth A. The cross-section 
of such a channel is also on the atomic scale, much less 
than £ and A. The theory of current transport through 
such short and narrow channels, often called quantum 
point contacts (QPCs), is well developed. For the single 
mode QPC, the Josephson current is given by 

eA Dsintj) , A J I - D sin 2 (0/2) 
I = ^^^^=^^^^= tanh — 

2h y/l-D sin 2 (0/2) fcsT 

(2) 

where D is the channel transparency and (f> is the phase 
difference. [39| 

The theory of nonstationary properties in SQPCs was 
developed by Averin and Bardas.[40] It is based on the 
idea of multiple Andreev reflections (MAR) [4l[ which re- 
veal themselves as current steps at subharmonics of the 
gap voltage 2 A/ en, at n = 2,3,.... They found that 
I-V characteristics strongly depend on the SQPC trans- 
parency and demonstrate the most pronounced subhar- 
monic steps at intermediate values of D ~ 0.4 to 0.5. The 
theoretical I-V curves reduce to a pure tunneling char- 
acteristics at D < 0.1 and to a featureless curve typical 
for an I-V of a clean S-N interface at D ~ l.[40j 



A7 SS - G NN A^ (3) 
2e 

where Gnn is the specific normal-state tunneling conduc- 
tance. 

It is easy to verify that the experimentally observed 
A7 SS in Fig. [2] is very close to the value given by Eq. 

The independence of A7 SS of electric stress indi- 
cates that the area of the tunnel barrier and its con- 
ductance does not change noticeably. It means that the 
total area of additional conduction channels formed by 
stressing is much smaller than the junction area. There- 
fore, a stressed junction can, in the first approximation, 
be represented by a parallel combination of the initial 
(unstressed) junction and some number of SQPCs (see 
Fig. El): I{V) = I (V) + 7 S qpc(V). To obtain the I- 
Vcharacteristics of these SQPCs, 7sqpc(^) we can sub- 
tract the I-V curve of the tunnel (unstressed) junction 
Iq(V) from the experimental I(V) in the stressed junc- 
tion. The result of this procedure is shown in Fig. [6j 
along with the theoretical I-V curves calculated for a 
SQPC with different transparencies, [iol] 

Although the subtraction procedure is not very accu- 
rate near the gap voltage because of the gap smearing 
in the highly stressed junction, it is clear that the I-V 
curve of the SQPCs has no significant current step at 
V = 2A/e. This is true for any stress current 1$ (any 
number of SQPCs) because in all the I-V curves in Fig. 
[2] the current step at V = 2A/e is the same as in the ini- 
tial junction. In the microscopic theory of nonstationary 
properties of a single SQPC,[4(| the size and broadening 
of the current step at V = 2A/e in the J-Vcurve as well 
as the size and broadening of subharmonic current steps 
strongly depends on the channel transparency. The step 
at 2A/e basically disappears at D > 0.6 while subhar- 
monic steps basically disappear at D > 0.9. Therefore, 
in the theory, the range of channel transparencies where 
the I-V curves would look like in Fig. [5£b) and Fig. O 
is very narrow 0.6 < D < 0.8. Hence, it would not make 
a significant error if, for simplicity, we assume that all 
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the channels have the same transparency. A better ap- 
proach would only be a full fitting of the I-V curves to 
the theory, which will be done elsewhere. It can be seen 
from Fig. [5] that all the I-V curves of SQPCs formed 
in the stress current range from the breakdown Tsb to 
2/sb lie between the theoretical curves with D = 0.6 and 
D = 0.8. Therefore, the average channel transparency 
formed at electric stressing can be taken as (D) = 0.7. 

From the linear part of the I-V curves at V > V e , 
the normal-state conductance of the channels is simply 



G 



SQPC 



R n0 , where i? n o is the normal-state 



tunnel resistance in the initial junction measured before 
stressing. Since the normal-state conductance of a single 
QPC is given by 2e 2 D/h, we can estimate the number 
N of created SQPCs using N ■ (D) = G SQ p C h/2e 2 . Its 
dependence on the stress voltage is shown in Fig. [9l At 
(D) = 0.7, the number of channels formed at the maxi- 
mum stress currents studied is ~ 2 x 10 3 , corresponding 
to an average channel density of 1.6 x 10 10 channels/cm 2 
and an average channel spacing of ~ 80 nm. In our model 
each channel is associated with a displaced ion in the bar- 
rier, so the estimates above give the surface density and 
the average distance between displaced ions in the bar- 
rier at the maximum stress currents = 100 mA. The 
obtained estimates are self-consistent with the initial sup- 
position that the number of additional conduction chan- 
nels is small and that the area occupied by the channel is 
negligible with respect to the area of the tunnel barrier. 

Since each SQPC carriers a supercurrent, the increase 
of the critical current in the stressed junctions A/ c = 

— IcO should be proportional to the number of conduct- 
ing channels. This is shown in Fig. [5] Using the assump- 
tion of equal channel transparency and (D) ~ 0.7, we can 
estimate the average contribution of each created channel 
SI C = AI C /N from Fig. Ej]to be about 0.115 ^A/channel. 
Using the experimental value of A/e = 1.35mV and 
(D) = 0.7, from Eq. we obtain the theoretical value 
of 5I C = 0.147 /zA/channel, a very close value. 

The theory also explains the existence of the excess 
current in the I-V curves as a nonzero average of the 
Josephson ac current at finite voltages. It can be seen in 
Fig. Othat the part of experimental I-V curves at V > V g 
(region of the excess current) agrees very well with the 
theory. The excess current defined as A/ cx = I cx — J ex o is 
shown in Fig. [5J where I ex was determined by fitting the 
I-V curves at V > V g to I(V) = I cx + V/R n and index 
identifies the initial, unstressed junction. The straight- 
line fit yields SI ex = 0.084 //A/channel, or 8I QX /I C = 
0.73 in very good agreement with the theory, [40j giving 
SI CX /I C = 2/tt at D ~ 1. 

So, by analyzing the additional subgap conductance, 
additional normal state conductance, and additional 
Josephson currents induced by electric stressing of the 
junction's barrier we obtained a self-consistent picture. 
In this picture, after a point of soft breakdown has 
been reached, atomic-size conduction channels with sim- 
ilar transparencies are being formed in the barrier as 
the stressing progresses. We also proposed a scenario 
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FIG. 7: Circuit diagram of the proposed model. In this model, 
the the application of electrical stress results in the formation 
of few additional conduction channels Gsqpc in the tunnel 
barrier which remains largely unchanged and is assumed to 
be the same as the initial, unstressed junction. 
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FIG. 8: The change in critical current and excess current 
caused by electric stressing as a function of the change in the 
normal-state conductance (bottom scale) and the number of 
created channels (top scale) based on (D) = 0.7. The straight 
line is the linear fit giving the average I c per channel of 0.115 
fiA and SI CX = 0.084 /iA/channel. 



in which the formation of these channels is a result of 
electric-field-stimulated ion migration away from the bar- 
rier. Let us see if this scenario also agrees with the avail- 
able experimental data. 

The treatment of ion transport in an electric field (elec- 
tromigration) is very well known and goes back to the 
original Frcnkel defect theory and works by Mott,[42j and 
Cabrera and Mott.[43[ Following [43|, the ionic current 
density can be written as 



Ji = 2anvq 



-(W~qaE)/k B T 



-(W+qaE)/k B T 



2anve~ w/ksT smh(qaE/k B T) 



(4) 



where n is the concentration of the mobile ions, v is the 
attempt frequency, W is the activation energy, q is the 
ion charge and a is the activation or half-jump distance. 
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FIG. 9: The number of created channels N (assuming (D) = 
0.7) versus the maximum applied stress voltage. The solid 
curve is the fit to Eq. (0) yielding Vo = 0.045 V. This is very 
close to the calculated value of Vo — 0.046 V. 

Ionic (atomic) transport in solids is usually proceeds via 
vacancy exchange mechanism (in which case W is the 
height of the energy barrier the ion needs to overcome) 
or via formation of Frenkel defects (vacancy-interstitial 
pairs) in which case W is the activation energy for forma- 
tion of a Frenkel defect. The total number of displaced 
ions in the barrier can be hence presented as 

N = At sinh(g AE/k B T) = At smh(V/V ) (5) 

where t is the stress duration, V is the voltage across the 
barrier, Vq = k^Td/qa, and d is the barrier thickness. 

Since in our model each displaced ion represents a new 
conduction channel, in Fig. [5] we plotted the "effective" 
number of conduction channels determined from the in- 
crease in the junction conductance N- (D)= GsQpc^/2e 2 
as a function of maximum voltage across the junction 
which develops during each stressing at room tempera- 
ture. The dependence given by Eq. is also shown and 
fits the data very well, giving V = 0.045±0.005 V. As- 
suming that oxygen ions are the mobile species, we can 
take q — 2e, a — 2.8A (nearest-neighbor distance), d = 1 
nm, and T = 300 K, and get from Eq. © V = 0.046 V in 
excellent agreement with the experimental data in Fig. [5] 
If instead we use the value qa = 4.28 eA determined from 
direct measurements of ionic current in alumina films, [37] 
the fit to the experimental data in Fig. [9] gives a tunnel 
barrier thickness d = 0.76 nm which is in a better agree- 
ment with the barrier thickness data (2f| [26| obtained 
from differential conductance vs voltage dependences. 

A secondary objective of this study was to relate 
the studied above effect of dc electrical stress with 
anomalous I-V characteristics of Josephson junction ar- 
rays frequently observed for many fabrication processes 
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FIG. 10: I-V curve of an as-fabricated 20-junction series ar- 
ray. Two junctions with significantly higher I c than the rest 
of the array are shown by arrows. Whereas the characteris- 
tics of these two junctions are qualitatively similar to those of 
the electrically-stressed junctions in this study, the amount of 
electric current that is required to damage the stressed junc- 
tions (7sb x(no. of junctions)) can not be supplied by the 
plasma processes employed during the fabrication. Hence, 
this fabrication-induced variation in I c in this array is un- 
likely to be caused by electron current-induced breakdown of 
the tunnel barriers. 



[20l l2ll . [44| . As was described in Sec. HI usually one or 
two junctions at the ends of the array demonstrate an I c 
that is significantly larger than the average I c of the other 
junctions in the array, see Fig. 1101 Clearly two junctions 
have an I c that is significantly larger than the average I c 
of the rest of the junctions. Usually, the junction with 
the highest deviation from the average is the last junc- 
tion in the array, the one which has the base electrode 
directly connected to the ground plane layer. [20L [2ll] The 
second deviating junction is the first junctions in the ar- 
ray, the one which is connected by the wiring layer to the 
chip contact pads. [13] These usually large deviations in 
I c are also reflected in the so-called knee region of the I-V 
curve where the knee current of the junctions in question 
is much higher and can be distinctly seen in the normal 
branch of the curve. Also the subgap conductance cor- 
responding to these two junctions is significantly higher 
than the average. MAR steps in the subgap region of 
these junctions have also been observed. [21| 

Qualitatively, these anomalous junctions look like the 
junctions subjected to electric stress as described above. 
Electric currents can flow through some of the Joseph- 
son junctions during their fabrication, e.g., during re- 
active ion etching. The question is whether the ampli- 
tude of the current is high enough to cause degradation 
of the barrier in many junctions on the wafer. From 
the results of dc electric stressing, we have found that 
the typical breakdown current in the smallest junctions 



10 



used in superconductor integrated circuits is 7 SB ~ 32 
mA at which a critical voltage develops across the bar- 
rier, i.e., Isb/Ic ~ 250. The typical 0.25-cm 2 test chip 
contains ~ 20 Josephson junction arrays. Therefore, to 
cause breakdown of all 20 grounded junctions in these 
arrays, it would need a current density of ~ 2.6 A/cm 2 
because for plasma-induced current the grounded junc- 
tions are connected in parallel. There are also several 
test chips spread out on the 150-mm wafers. So the total 
current through the wafer should be huge. The typical 
rf power used for etching in our process is from 40 to 
150 W and peak-to-peak voltage is from 305 V to 900 V 
depending on the circuit layer, corresponding to power 
density from 88 mW/cm 2 to 332 mW/cm 2 . Although, 
rf power is coupled capacitively to the wafer and to the 
circuits' ground plane, the rf current can flow directly 
through the junctions which are galvanically coupled to 
the ground plane. Estimating the amplitude of electric 
current supplied by the rf source, we find it to be in the 
range of 0.5 A to 0.7 A, corresponding to the current den- 
sity from 1 mA/cm 2 to 1.55 mA/cm 2 . This is clearly too 
small a current to cause significant damage to oxide barri- 
ers in Josephson junctions. Therefore, we conclude that, 
although exposure to processing plasma has a potential 
in creating voltage differences on the wafer surface which 
are large in comparison with the typical breakdown volt- 
ages in ultrathin oxide barriers, [45j the typical etching 
plasmas used in superconductor integrated circuits fabri- 
cation do not supply enough current to electrically break 
down the tunnel barriers of Josephson junctions due to 
their high tunneling conductance. The same conclusion 
was also reached in 



21] 



Because of these results, it is likely that the damage 
to junctions galvanically coupled to the ground plane or 
other circuit layers observed in [2(| is chemical or elec- 
trochemical in nature. For instance, it can be due to 
diffusion or electromigration of impurity atoms from the 
layers in direct contact with the junction to the junction's 
electrodes and to the barrier, e.g., hydrogen atoms dis- 
solved in Nb. In many respects the features of oxide bar- 
rier damage by H + electromigration should be similar to 
that of electron current damage. However, the required 
currents can be substantially lower because hydrogen is 
known to be highly mobile in Nb. Therefore, it should 
be much easier to cause diffusion of hydrogen atoms in 
Nb than to displace oxygen atoms from the AlO^ ox- 
ide barrier. This damage mechanism will be considered 
separately. 



VI. CONCLUSION 

We have studied the effect of dc electric stress (current) 
applied at room temperature to Nb/Al/AlO^/Nb junc- 
tions with J c = 1.0 and 4.5 kA/cm 2 on their Josephson 
and quasiparticle tunneling properties. We have observed 
a very soft breakdown: above a certain threshold stress 
current Isb the changes in I-V characteristics indicate 



formation in the barrier of additional conduction chan- 
nels with increased transparency. These channels reveal 
themselves in dramatic increase in the subgap conduc- 
tion, appearance of subharmonic current steps (multiple 
Andreev reflections), an increase in the critical current, a 
decrease in the normal state resistance, and an increase 
in the I c R n product. The range of the stress currents 
where these gradual changes occur is quite broad, from 
-^SB to > 2/gB, corresponding to ~ 250J C to > 500/ c . The 
threshold (breakdown) current depends on the stress cur- 
rent polarity and is larger when the Nb counter electrode 
is at a positive potential. 

We have suggested that the soft breakdown proceeds 
via formation of defects in the barrier as a result of ion 
migration in the applied electric field. The stress polar- 
ity dependence was explained as due to internal electric 
field caused by the difference in the work functions of 
the junction electrodes. The observed breakdown fields 
E c ~ 3.85 x 10 8 V/m {V c S 0.385 V across the barrier) 
correspond to the magnitude required for near-neighbor 
displacement of oxygen ions. 

From the changes in the tunneling characteristics in the 
superconducting state we concluded that the total area of 
the additional conduction channels is much smaller than 
the area of the initial tunnel barrier, and that they can be 
considered as connected in parallel to the initial tunnel 
junction whose properties remain virtually unchanged by 
electric stressing. Using this parallel connection model, 
we estimated the transparency of the additional conduc- 
tion channels by comparing their I-V characteristics to 
the nonstationary theory of current transport in super- 
conducting quantum point contacts by Averin and Bar- 
das. We found that all the channels have similar trans- 
parencies D in the range from ~ 0.6 to ~ 0.8. We found 
an excellent agreement of the shape of I-V characteristics 
of the conduction channels created by electric stressing 
with the Averin-Bardas theory for SQPCs.[40[ From the 
change in the normal-state resistance we calculated the 
effective number of the additional conduction channels 
N. The observed changes in the junction critical cur- 
rent per channel and the excess current per channel were 
found to be consistent with the predictions of the micro- 
scopic theory for SQPCs. 

We also found an exponential dependence of the num- 
ber of additional conduction channels created by electric 
stressing on the stress voltage, N ~ sinh(V/Vb). The 
characteristic voltage Vq was found to be in excellent 
agreement with the proposed model of defect formation 
by ion electromigration out of the barrier. 

Due to the high tunneling conductance of oxide barri- 
ers used in superconductor digital electronics it is unlikely 
that breakdown currents can be reached during plasma 
processing steps of junction fabrication. However, it is 
possible that electromigration and diffusion of impurity 
atoms in Nb layers is responsible for pattern-dependent 
tunnel barrier degradation frequently observed in super- 
conductor integrated circuits. 
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